Finite temperature effects on the collapse of trapped Bose-Fermi mixtures 
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By using the self-consistent Hartree-Fock-Bogoliubov-Popov theory, we present a detailed study 
of the mean-field stability of spherically trapped Bose-Fermi mixtures at finite temperature. We 
find that, by increasing the temperature, the critical particle number of bosons (or fermions) and 
the critical attractive Bose-Fermi scattering length increase, leading to a significant stabilization of 
the mixture. 
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The recent experimental achievement in the trapping 
and cooling of degenerate Bose-Fermi mixtures of alkali- 
metal atoms have introduced interesting new instances 
of a quantum many-body system 0, 0, 01 ■ Similar 
to the purely Bose gases such mixtures provide a unique 
opportunity of investigating fundamental quantum phe- 
nomena. Up to date a number of theoretical analyses of 
trapped Bose-Fermi mixtures have been presented, ad- 
dressing, for example, the static property Q, the phase 
diagram M], the behavior of low- lying collective excita- 
tions H H EH [H [H [H El, the expansion 0, as 
well as the relevant implications for the achievement of 
Bardeen-Cooper-Schrieffer (BCS) transition to a super- 
fluid phase of the fermionic components 0, 0, . The 
last point is particularly interesting and is highly favor- 
able by tuning the Bose-Fermi interaction by means of 
Feshbach resonance |2(j , in order to induce a large 
effective attraction between fermionic components by ex- 
changing density fluctuations of the bosonic background 
\wL Il7l llSj . It has been suggested in Ref . [2(| that the 
optimal conditions for boson-induced Cooper pairing are 
reached when the Bose-Fermi mixtures are close to a crit- 
ical point, beyond which the system will become unstable 
and collapse. 

This kind of collapse is indeed observed very recently in 
a binary 40 K— 87 Rb mixture |2l|. The observation indi- 
cates a large attractive interaction between bosonic and 
fermionic atoms, and therefore highlights the possibil- 
ity of inducing a boson-mediated effective attraction be- 
tween two fermionic components once another hypcrfinc 
state of fermionic atoms is populated. On the theoretical 
side, the collapse or the instability of a trapped Bose- 
Fermi mixture at zero temperature has already been 
considered by several authors d IH |H HI |2g. In 
Refs. |H H, the mean-field instability of a spheri- 
cally trapped Bose-Fermi mixture is studied using a stan- 
dard two-fluid model for the condensate and degenerate 
Fermi gas. Most recently, the same mean-field model has 
been used in Ref. [25|. in which the precise geometry of 
the experiment is taken into account. Its prediction has 



been compared with the experimental stability diagram 
of 40 K- 87 Rb mixtures to give a better estimate of the s 
-wave Bose-Fermi scattering length 25] . 

In the present paper we investigate the effects of the 
temperature on the stability of a spherically trapped 
Bose-Fermi mixture by only a mean-field approach. The 
motivation is twofold: At first, in the experiment there is 
a detection limit for determining the temperature. This 
limit is around 0.67T C , where T c is the transition temper- 
ature of the Bose-Einstein condensation (BEC). Below 
this characteristic temperature, though no uncondensed 
fraction of the BEC is practically detectable, the possible 
role played by the thermal cloud of bosons in the stability 
of the mixture should be clarified. On the other hand, 
in the context of boson-induced BCS superfluidity, the 
critical temperature for the formation of s-wave Cooper 
pairing increases exponentially as the mixtures move to 
the collapse point. The qualitative or quantitative under- 
standing of the collapse of Bose-Fermi mixtures at finite 
temperature is therefore crucial to optimize the experi- 
mental conditions to achieve the BCS transition. 

Our investigation is based on the simplest self- 
consistent mean-field theory — the Popov version of the 
Hartree-Fock-Bogoliubov (HFB) theory — that has been 
generalized by us to the system of Bose-Fermi mixtures 
in previous works |l4 . |26| . For a weakly interacting Bose 
gas, it was shown that the HFB-Popov theory gives the 
correct thermodynamic properties with a very good accu- 
racy j2]j , though it fails to predict the correct behavior of 
collective excitations at high temperatures [2^]. In Ref. 
|29) . this theory has been applied to study the collapse 
of attractive Bose-Einstein condensates. 

In the following we briefly summarize the main points 
of HFB-Popov theory. The trapped binary mixture is 
considered as a thermodynamic equilibrium system un- 
der the grand canonical ensemble whose thermodynamic 
variables are Nb and Nf, respectively, the total number 
of trapped bosonic and fermionic atoms, T, the absolute 
temperature, and /i& and ///, the chemical potentials. In 
the second quantization language, the density Hamilto- 
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nian of the mixture reads(?i = 1), 



where the operator 
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where f/K 1 ") (</K r )) is the Bose (Fermi) field operator that 
annihilates an atom at position r. The first and sec- 
ond terms in the brackets contain the kinetic-energy op- 
erators and the external trapping potentials V^ xt (r) — 
mbUJ b 2 r' 2 /2 and V/ Xt (r) = nifwir 2 /2 for the bosonic 
and the fermionic species with masses nib and m/, re- 
spectively. In the dilute regime, we have considered 
two types of contact interactions: the interactions be- 
tween bosons, and the interactions between bosons and 
fermions. They are parametrized by the coupling con- 
stants gbb — Airfi 2 abb / rrib and gbf = 2i:h 2 abf /m r , respec- 
tively, in terms of the s-wave scattering length abb and 
a b f, withm r = m b m f I (mb+m /) being the reduced mass. 
We have neglected here the fermion-fermion interactions 
since for spin-polarized fermions the s-wave contact in- 
teraction is prohibited by the Pauli principle. The next 
leading order, p-wave interaction is small at low energy 
[3fj| and will not be considered in the following. 

We consider the many-body ground state | Vf) of Hamil- 
tonian (JTJ as a direct product of a symmetric iVf,-body 
state \^b) for the bosonic species and an antisymmetric 
Nf -body state \^ f) for the fermions. With this choice we 
are not considering the possible correlation effects beyond 
the mean-field approximation [3lj . The density Hamilto- 
nian describing the Bose-Fermi coupling can therefore be 
decoupled in a self-consistent mean-field manner, namely, 

Hit - 9bf + (^)^- ^ + ^){<P + 4>)} ■ 

(2) 

As a result, the bosonic and fermionic atoms experience, 
respectively, the effective potentials V^Ar) = V^ xt (r) + 

5b/n/(r) andV e ^ / (r) = V e f xt (r)+g b fn b (r) 7 with n b (r) and 
n/(r) being the bosonic and fermionic density distribu- 
tions, respectively. It is then straightforward to apply the 
HFB-Popov theory for the bosonic species, following Ref. 
32]. After invoking the Bose symmetry breaking in the 
equation of motion for the Bose field operator ip(r), we 
obtain, respectively, the modified Gross-Pitaevskii (GP) 
equation for the condensate wave function <&(r), 



£ GP $(r) = 0, 



(3) 



and the modified Bogoliubov-deGennes (BdG) equations 
for the thermal quasiparticle amplitudes Ui(r) and i>i(r), 

[C GP + g b bn c (r)} itj(r) + g 6b n c (r)w 4 (r) = ejWj(r), 
[£ G p + gbbnc{r)]vi(r) + g bb n c (r)ui(r) = -e;Ui(r),(4) 



GP 



V 

2rrib 



- VLt - f^b+9bb (n c (r) + 2n T (r) ) +g bf n f (r) , 

(5) 

is the generalized GP Hamiltonian that includes the 
mean-field contributions generated by the interaction 
with the thermal cloud and the fermionic species. Once 
these wave functions have been determined, the local den- 
sity of the condensate and of the noncondensate, and the 
total bosonic density distribution can be calculated ac- 
cording to 



n c (r) 
n T (r) 
n b (r) 



|$(r) 



E 



Ui (r)| 2 + Mr)| 2 



e /3e 4 

n c (r) +n T (r), 



(0) 



with (3 = 1/ksT being the inverse temperature. To solve 
the modified GP and BdG equations, one has to evalu- 
ate the fermionic density distribution nf(r). To this aim, 
we employ the finite-temperature Thomas-Fermi approx- 
imation (TFA) pj . 
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Specially, for the range of Nf considered here (Nf ~ 10 4 ), 
it was shown that TFA is a good approximation for a mix- 
ture with repulsive Bose-Fermi interaction at zero tem- 
perature 0- In our case, this approximation is expected 
to work even better, since the quantum fluctuations are 
suppressed by the finite temperature and the increased 
density due to the Bose-Fermi attraction. 

Equations ©-© form a closed set of equations that 
we have referred to as the "HFB-Popov" equations for 
a dilute Bose-Fermi mixture. The simultaneous solution 
of this set of equations gives the temperature-dependent 
density profiles of the condensate, of the noncondensate, 
and of the degenerate Fermi gas, from which we can ex- 
tract the stability conditions of the system. We have nu- 
merically solved these coupled equations by an iterative 
scheme, as described in Ref. [26j. 

With above tools we investigate the effects of the tem- 
perature upon the collapse of a trapped Bose-Fermi mix- 
ture induced by the attractive Bose-Fermi interaction. 
Rather than taking the realistic experimental conditions, 
in this paper we present a general analysis of finite tem- 
perature effects, and restrict ourselves to spherical sym- 
metric systems. The reason is that for this geometry the 
temperature-dependent density profiles of the system are 
easily calculated. On the contrary, for anisotropic sys- 
tems, though the formalism given in the present paper is 
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FIG. 1: Radial density profiles of a Bose- Fermi mixture with JVj, = Nf = 5 x 10 4 for different temperatures and Bose- Fermi 
interaction strengths. The condensate and noncondensate densities n c (r) and nr{r) are given, respectively, by the solid and 
dash-dotted lines (left scale), and the fermion density nf(r) by dashed line (right scale). To visualize riT{r) is enlarged by a 
factor of ten. The upper and lower rows show, respectively, examples with decreasing Bose-Fermi attraction at T /Tc = 0.005 
and 0.67, where Tc = 0.94?kjiV 6 1/ ' 3 /ks is the transition temperature of a non-interacting Bose-Einstein condensate in the 
thermodynamic limit. We fix the boson-boson scattering length abb /I = 0.005. 



suited as well, the numerical calculation is computation- 
ally much more involved. In particular, for the highly 
anisotropic traps used in the experiments at the Eu- 
ropean Laboratory for Non-linear Spectroscopy (LENS) 
0, 0] it is unlikely to obtain useful information without 
the use of further approximation. In the following, we 
consider equal masses for the two species m = = m,f 
and identical trapping frequencies uj = u>b = u>f. The 

1/2 

oscillator length I = (h/mui) and trap energy hu>, 
respectively, serve as fundamental length unit and en- 
ergy unit. We also take the transition temperature of 
a non-interacting Bose gas with Nb = 5 x 10 4 atoms, 

1/3 

as a character- 
istic temperature. In most calculations, we have set 
T = 0.67T C that corresponds to a typical detection limit 
for the measurement of the temperature. 

First we consider the effects of the temperature upon 
the density profiles of the mixture. The upper and lower 
rows in Fig. 1 show the density profiles of configura- 
tions with Nb — Nf — 5 x 10 4 particles for three dif- 
ferent values of the s-wave Bose-Fermi scattering length 
a bf /l at two temperatures: T = 0.005T C and T = 0.67T C . 
The value of abb /I = 0.005 has been fixed, which corre- 
sponds to abb ~ lOOaBohr for a typical trap with I = lfim. 
Our results at the lower temperature T = 0.005T C are 
in good agreement with the findings by Roth and co- 
workers [23, 0| > that is> the mutual attraction between 
bosons and fermions results in an enhancement of both 
densities in the overlap region. In particular, as shown in 
Fig. (lc), both the bosonic and fermionic densities grow 



substantially when the mixture is close to the instabil- 
ity point (note that at T — 0.005T C the critical value of 
abf/l is —0.0193). This enhancement, however, is much 
reduced at a finite temperature, whose effect is a broad- 
ening of the density distributions of the condensate and 
of the Fermi gas and therefore reduces their center den- 
sities. As can be seen in Fig. (If), at T = 0.67T C both 
densities is decreased by a factor of 2 compared to the 
lower temperature case. As a consequence, the mixture 
in this case is expected to be much stabilized against 
collapse. 

In a simplified model , the collapse or the instability 
of the mixture is governed by the balance between the 
kinetic energy of fermions and the mutually attractive 
mean-field generated by the Bose-Fermi interaction (see, 
for example, the discussion above the Eq. (11) in Ref. 
6J), If the Bose-Fermi attraction becomes too strong, 
i.e., the numbers of bosons (or fermions) or the scatter- 
ing length abf becomes too large, the attractive mean 
field cannot be stabilized by the kinetic energy anymore, 
so both density distributions grow indefinitely within the 
overlap region and collapse. According to the study re- 
ported in Refs. |(| |23, yj, |2||, the onset of this mean- 
field instability is monitored by the failure of the conver- 
gence during the iterative procedure. In this manner, we 
can determine the critical s-wave Bose-Fermi scattering 
length abf or critical particle numbers jV^™* and Nf nt 
beyond which the collapse occurs. 

In Fig. 2, we show the critical particle number as 
a function of the s-wave Bose-Fermi scattering length 
abf/l with fixed abb/l = 0.005 for two temperatures. In 
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FIG. 2: Behavior of the critical particle number as a func- 
tion of the s-wave Bose- Fermi scattering length atf/l for 
N b = Nf at two temperatures: T/T c = 0.005 (solid line) 

and T/Tc = 0.67 (dashed line). Tc = 0.9ihujN b 1/3 /k B is the 
ideal transition temperature evaluated at Nb = 5 x 10 4 . The 
boson-boson scattering length abb /I = 0.005 is fixed. 
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FIG. 3: (a) The critical number of bosons N§ r% as a func- 
tion of the rescaled temperature T/Tc with fixed number of 
fermions Nf = 5 x 10 4 . The dashed line shows the correspond- 
ing number of condensed bosons, (b) The critical number 
of fermions Nf™* as a function of the re-scaled temperature 
T/Tc with fixed number of bosons Nb — 5 x 10 4 . In both cases 
Tc = 0.94hu!N^ 3 /ks is the ideal transition temperature for 
Nb = 5 x 10 4 . The other parameters are: abb/l = 0.005 and 
dbf/l = -0.020. 
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FIG. 4: Region of stability of the Bose-Fermi mixture, as a 
function of the number of atoms, for the case of abb/l = 0.005 
and abf/l = —0.020. Lines show the boundary between the 
stable (left) and collapse (right) regions, for three values of 
the temperature: T/T c = 0.005 (solid line), T/T c = 0.50 
(dashed line) and T/Tc = 0.67 (dash-dotted line), where 
Tc = 0.94:hujN^ 3 /Ub is the transition temperature for a non- 
interacting Bose gas with Nb — 5 x 10 4 . 



this calculation we have kept equal numbers of bosons 
and fermions (N b — Nf). At the lower temperature 
T = 0.005T C , we observed that even a small decrease 
of Bose-Fermi attraction can reduce the critical particle 
number significantly, i.e., N cnt changes by a factor of 2 
when abf/l changes by 5 percent. This strong dependence 
of the criticality on abf/l is in accordance with the scaling 
law for the critical particle number: N crlt ~ |afc/|~ a with 
a = 12, as discussed in Refs. [E E2> EH • The presence of 
a moderate temperature T = 0.67T C results in a substan- 
tial increase of the critical particle number. For example, 
at a bf /l = -0.022, N crlt grows from 1.1 x 10 4 to 3.6 x 10 4 
with the inclusion of the temperature. Apart from this 
growth, the critical particle number is still strongly de- 
pendent on abf/l. However, the presence of the finite 
temperature leads to a re-normalization of the scaling 
exponent, i.e., a ~ 8. 

In Figs. (3a) and (3b), we study the critical par- 
ticle numbers of bosons and fermions as a function of 
the temperature. We report, respectively, the predic- 
tion on the critical particle numbers of bosons (with fixed 
Nf — 5x 10 4 ) and fermions (with fixed Nb = 5x 10 4 ). For 
comparison, in Fig. (3a), we also show the critical num- 
ber of the condensed atoms against temperature(dashed 
line). As expected, the critical particle number for each 
species increase with increasing the temperature. The 
dependence is highly nonlinear. Below 0.5T C , the criti- 
cal particle number varies slowly with the temperature, 
whereas above 0.5T C it rises up steeply. In addition, the 
critical particle number of fermions grows more rapidly 
than that of bosons against the temperature. 

Finally, we have built a region of stability of the Bose- 
Fermi mixture in Fig. 4, as a function of the number of 
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atoms, for the cases of abb/l — 0.005 and a^f/l — —0.020, 
for three values of the temperature: T/Tq = 0.005, 
T/T c = 0.50 and T/T c = 0.67. Each of the curves marks 
the limit of stability. For numbers of bosons or fermions 
below the stability limit the mixture is stable, otherwise 
the mixture is unstable against mean-field collapse. Fig- 
ure 4 indicates that the region of the stability broadens 
with increasing temperature. This behavior emphasize 
again that the inclusion of a temperature gives rise to a 
significant stabilization of the mixtures. 

In conclusion, we have investigated the mean-field sta- 
bility of a spherically trapped binary Bose- Fermi mixture 
at finite temperature. We solved the coupled HFB-Popov 
equations numerically and obtained the critical particle 



number as a function of the temperature and the s-wave 
Bose-Fermi scattering length. We have shown that the 
critical particle number and the critical Bose-Fermi scat- 
tering length increase with the inclusion of a moderate 
temperature that corresponds to the typical experimental 
detection limit. This leads to a significant stabilization 
of the Bosc-Fcrmi mixtures. 
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